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Using the valence-bond-solid (VBS) approach and the Schwinger boson mean field 
approximation, we study the dependence of the Haldane gap of a spin-1 linear chain 
Heisenberg antiferromagnet on impurity doping with different spins. The impurity 
spins affect the singlet pairing order parameter A and the constraint factor A. As 

2/3 

a result, the Haldane gap is reduced by a factor ~ , with n« as the impurity 
concentration, and eventually collapses at Hi ~ 1/t; with £ as the VBS correlation 
length. This theoretical prediction can be verified by neutron scattering experiments. 
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Some years ago Haldane JTJ conjectured that the excitation spectrum of a linear-chain 
Heisenberg antiferromagnet (LCHA) with integer spin has a finite energy gap Eh above 
its singlet ground state, while a LCHA with half-integer spin has a gapless spectrum. This 
conjecture has been strongly supported by an exact solution of a specific model 0, numerical 
studies 0, and experiments 0,0. Recently, there has been a revived interest in this 
problem due to the experimental observation of spin 1/2 degrees of freedom 0, at the 
ends of a finite spin-1 chain induced by doping with different spins, as predicted by theory 



P.W. Anderson [IU| has proposed the Resonant- Valence- Bond (RVB) model to describe 
some spin 1/2 antiferromagnet systems. Affleck, Kennedy, Lieb and Tasaki (AKLT) have 
generalized this idea to spin-1 LCHA by proposing the Valence-Bond- Solid (VBS) state. In 
this VBS state, the valence bonds formed by two 1/2 spins as a singlet fj, — ||? connect two 
nearest neighbors, while the two 1/2 spins on the same site should be symmetrized to form 
a triplet state 5=1. This is a translationally invariant, singlet state. The spin correlation 
functions decay exponentially and there is a gap Eh in the energy spectrum which is closely 
related to the VBS order parameter. 

Consider a generalized Heisenberg Hamiltonian 

Hq = Jo^2^Si ■ S i+ i — f3(Si ■ S i+ i) 2 , (1) 

where — l</3<lisa parameter. The VBS state proposed by AKLT is an exact, non- 
degenerate ground state of ([!]) for (3 = —1/3. The theoretical analysis using the non- 
Abelian bosonization technique and the conformal field theory [PH| , as well as the exact 
diagonalization for small clusters seems to show that systems described by (|T|) with 
— 1 < (3 < 1 belong to the same universality class as the /3 = — 1/3 case, so its properties 



are generic to all systems of this class. Recently, Ng [fy| has applied the Schwinger boson 



mean field approximation (SBMFA) adopted by Arovas and Auerbach |T3| for quantum spin 



discussed in his paper and will not be repeated here. 

In this Letter we use the SBMFA to consider the effect of random doping with different 
spins on the Haldane gap. In this approach the Haldane gap Ejj is determined by the VBS 
order parameter A and the single-occupancy constraint factor A. Upon doping, apart from 
the end states, A and A are modified by both spin-flip scattering and the difference in the 
spin values. By solving self-consistently the system of equations for the order parameter and 
the constraint we find that the reduction of the Haldane gap Eh originating from the VBS 
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structure is proportional to n/ , where is the impurity concentration, and it collapses 
when rij is comparable to the inverse of the VBS correlation length. The situation here 
is rather similar to superconductors doped with paramagnetic impurities which break the 
time-reversal symmetry of the singlet pairing, giving rise to reduction and eventual collapsing 



of the energy gap ||14j| . Of course, if the spin chain is strictly one-dimensional, doping by 
non-magmetic impurities will break the chain. However, the real system is only quasi-one- 
dimensional, so the interchain coupling as well as the superexchan The dynamic structure 
factor S(q,u), measured in the neutron scattering experiments, is the Fourier transform of 
the spin-spin correlation function < Si(t) ■ Sj(0) >, where < • • • > means thermal average. 
Apart f 



The Hamiltonian of a LCHA doped with impurity spins can be written as fll6 
H — Hq + H im , H = J Yji Si ■ Si+i, H lTn = J2<ia> A a ■ Si, 

(2) 

A a = gStT ~ JoS a , S 2 = S(S +1), J > 0, 

where S im is the impurity spin, g is the coupling between impurity and host spins, a denotes 
impurity site, < ia > means summation over all impurity spins with % = a ± 1. 

We will discuss a generic case when impurities are randomly distributed along the LCHA 



and the averages are S im = 0, S 2im = 1/3 S im (S im + 1) (where 9 means any direction), 
A 2 = g 2 S im {S im + 1) + J 2 S(S + 1). 



[a, b'] = 0, [a,a^] = 1. It is convenient to use Nambu's four- component formalism to treat 
the Hamiltonian @. Let ipi = (a i? b u bj, aj) T , ip{k) = (a k , b k , b_ k , a)_ k ) T , then ipj = 
(af, b\, bi, di), ip\k) = (a\, b\, a_fc), where T means transposition and k is the wave 
vector. To carry out the decoupling in the SBMFA we introduce the VBS order parameter 
Aij =< ciibj — bidj >, with ij denoting the nearest neighbor pair of sites. 

We first consider the undoped system. In the SBMFA the Hamiltonian can be written 



as 



i i 

+ £J|A M+ i| 2 + £A i(-l-2S) 



(3) 



where 
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while * means complex-conjugate, J = Jo/2, / and 03 are 2x2 unit and Pauli matrix, 
respectively. 

Performing the gauge transformation of Read and Newns ]T7[ to absorb the phase factor 
and to remove the redundancy of constraint, the Lagrangian corresponding to eq.(^) can be 
written in the momentum space as 



L = 1/2^^^)^! - X Q 2 + JA u k Q 3 )ifj(k) + const, 



(4) 



where v k = 2 sin k, A is the module of the VBS order parameter. Here we take the lattice 
spacing a as the length unit and denote Q = (Qi, f2 2 , f^) as 
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From eq.(^j) one obtains the boson Green function 

uQi + A ^2 + JA i/ k Q s 



n.iu , a 



where 77 is a positive infinitesimal. The Schwinger boson excitations form a continuum band 
and their dispersion relation is e(k) = (A 2 , — 4J 2 Aq sin 2 k) 1 ! 2 . 

Now we can easily calculate the dynamic structure factor at T = OK as follows 



S(q,u) 



12tt 



dki 
2^ 



duo 2 



5(uJ — LO\ — UJ 2 ) 



(6) 



{ImG 044 (fci, u;i)ImGoii(g + h, co 2 ) - ImG 042 (-A; 1 , Ui)lmG 013 (q + k u u 2 )}, 

where G 0i4: ,G ii,G 0i2 and G013 are elements of the Green function matrix (||). Integrating 
over uji and uj 2 of eq.(^), we find that the 5-function representing the energy conservation 
becomes 8{u - (A 2 , - 4J 2 A 2 I sin 2 ki) l/2 - (Aq - 4J 2 A 2 1 sin 2 (g + ki)) 1/2 ). Hence only when 
^ > 2(Aq — 4J 2 Aq) 1 / 2 , a neutron can be scattered and there must be two Schwinger bosons 
to form a physical magnetic excitation for spin S = 1. It is also clear from the Schwinger- 
Boson constraint for the spin-1 chain (a^a + Wb = 2) that they must appear in pairs. 
This is entirely consistent with the description based on field-theoretical mappings [JTQj . 



Thus the Haldane gap Eh is the bottom of the spin excitation continuum and is given by 
Eh = 2(Aq — 4J 2 A 2 >) 1 / 2 . The parameters A and Ao can be consistently evaluated ]12|, ||IBf . 
At temperature T = OK, they are A = 4.7296J, A = 2.36312, E Ho = 0.1772J . 

Now we consider the doped system. In the following A a is referred to as a quasiclassical 
quantity. In the Nambu representation the impurity part of the Hamiltonian (|2|) can be 
written as H im = 1/4 ^ ipj (A a ■ E)^, in which E = (E 1 , E 2 , E 3 ) and 
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with oi, cr 2 , cr 3 as Pauli matrices. 

Since the impurities are randomly distributed, their average effect is zero if each impurity 
scatters a Schwinger boson only once. Thus the multiple-scattering effect should be consid- 
ered. Analysing the Feynman graphs, one finds that the contributions from the crossing dia- 



function G(k, u) of the doped system after averaging over the random distribution recovers 
the translational invariance. Hence we have u) = Gq 1 (k, to) — u), where S(A;, u) 

is the self-energy term due to the impurity scattering and Gq 1 (k, to) = luQi — WI2 + JAu^3, 
which has the same form as eq.([5]). However, parameters A and A should be self-consistently 
determined like Ao and A and in principle they are different from the latter due to the 
impurity effect. Furthermore, we make an ansatz that the inverse of the Green function 
G (k,z) = ZzQ\ — XQ2 + JAZi/fcf^, where z is the frequency in the complex plane, Z is 
the wave function renormalization factor, A and A are the renormalized parameters for the 
A-multiplier and the order parameter, respectively. All of them are functions of k and z. 
The self-energy part after averaging over random distribution can be written as 

(A ■ Z) 2 G (k u u) cos 2 (k-h)-±. (7) 

From the definitions given above and eq.(0), we obtain a closed set of integral equations to 
determine Z, A, and A. Since the doping density n« is rather small, they can be expanded 
in powers of rij and we will keep only the leading order of in the following. Accordingly 



we obtain that 1 18 



Zz = z(l -r{k,z)), A = A(l + r(k, z)), AZ = A(l + r(k, z)), (8) 

where 

r(k, z) = ^ I (4 cos 2 k + ~jP A l° S2k )(z 2 - \ 2 )~ 1/2 (z 2 - A 2 + 4J 2 A 2 ) _1/2 - ^ J . 

In this case A and A should also be evaluated consistently. If the impurity density rij is zero, 
one recovers Ao and Ao- When rii is very small, they can be expanded up to n.j. 

Now we consider the density of the Schwinger boson states for the doped system given 

by 

dN 1 /• T .dk 
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(10) 



and dN/dz is zero for c<j < oo min . Here = 2(A 2 — 4J 2 A 2 ) 1 / 2 . 

In SBMFA, the dynamic structure factor S(q,u) is a two-particle Schwinger- boson 
Green function. The random doping effect consists of two parts. The first is included in the 
simple product of two single-particle Schwinger-boson Green functions G(k,u). The second 
is the correlation effect between two Schwinger bosons induced by scattering at the same 
impurity. This correlation effect is mostly contributed by the two-particle Green function's 
ladder diagrams and maximum crossing diagrams. Their contribution is at least proportional 
to the impurity density n{. In the weak doping case we may ignore this effect. Within this 
approximation S(q,u) in the presence of impurity scattering can be easily calculated by 
replacing G in eq. (|5p with G renormalized by impurity scattering. As in the undoped case, 
S(q, uj) vanishes for is the renormalized Haldane gap of the doped sample. If we keep eq. (pTo|) 
up to the first order of nj, Eh and A can be replaced by Eh and Ao, respectively. Now we 
make a numerical estimate for the case of non-magnetic doping. Such an experiment has been 
performed by Glarum et al. || For such non-magnetic case we have A 2 = JqS(S + 1) = 8J 2 
(for S = 1 ). Substituting the parameter values for A and A into (0), we find that 



from which we see that En ph when rii ph 7.8%. In the SBMFA treatment the VBS 



correlation length is comparable with the average impurity distance. This is reasonable and 
could be expected from our argument at the very beginning. 

The doping dependence of the Haldane gap can be observed directly in neutron scattering 
experiments, for example, in NENP doped with Cu. It will also show up in other physical 



E H = E Ho (l - 5.6722n- /3 + 0.4775n0 



(11) 



state's correlation length £ is about 10 lattice spacings [T2|], [jC^ , i.e., almost twice the exact 
value Q. Within the same approximation we find that the Haldane gap collapses when the 



available up to now. However, it has been observed that the spin 1/2 ESR signal decays very 
rapidly at much lower temperatures than the Haldane gap itself 0, ||. Mitra, Halperin and 
Affleck |20f have suggested to explain this effect by considering the thermal excitations which 
will shift the ESR frequency out of the observation window for short and intermediate chains, 
while broadening the signal for very long chains. This effect should be there. However, we 
would suggest that there is a significant reduction of the Haldane gap itself (about 30 %) 
even at 1 % doping by non-magnetic impurities due to the effect considered in this Letter. It 
seems to us that these two complementary aspects should be both taken into account while 
interpreting the data. Our calculation was done for zero temperature, but it can be easily 
extended to finite temperatures to compare with actual experiments [|T^]. It is known |2T| 
that the Haldane gap for undoped systems first incr 

In conclusion we have considered the effects of doping by impurity spins in LCHA on 
the Haldane gap. In this calculation we have considered only the impurity effect on the 
Schwinger boson continuum spectrum which in turn leads to reduction of the Haldane gap 
itself, like the early calculation of Gor'kov and Abrikosov [111] for the paramagnetic impurities 



in superconductors. However, there are bound states in both spin chain problem ||, |T2 
and superconductors [Q. The bound states in superconductors are rather shallow and it is 
difficult to observe them experimentally. On the contrary, the bound states in the LCHA 
are rather deep, very close to the ground state. Hybridizing with the ground state, they give 
rise to the spin 1/2 signal observed in experiments. They are localized within the correlation 
length and in the low low doping case we can still focus on the scattering states. On the 
other hand, if we discuss the collapse of the Haldane gap upon doping quantitatively, both 
tail states of the continuum band and the band formed by these localized states should be 
taken into account. In this case the impurity density rij is not very small any longer and we 
cannot neglect the correlation effect between the Schwinger bosons coming from scattering 
on the same impurity. In fact, the properties of such gapless LCHA with nonvanishing 



impurity in the normal metal and superconductors. Their quantum fluctuations should have 
important effects. The work on these issues is in progress and will be reported later. 
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